Orbitally Degenerate Spin-1 Model for Insulating V2O3 
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Motivated by recent neutron, X-ray absorption and resonant scattering experiments, we revisit the 
electronic structure of V2O3. We propose a model in which S=l V 3+ ions are coupled in the vertical 
V-V pairs forming two-fold orbitally degenerate configurations with S=2. Ferro-orbital ordering of 
the V-V pairs gives a description which is consistent with all experiments in the antiferromagnetic 
insulating phase. 
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Although the metal- insulator transition in V2O3 has 
long been studied as a classic Mott-transition [1-3], the 
detailed electronic structure remains open. Recently 
new experimental techniques have been applied but these 
have not resolved the issue. Rather they have reopened 
the long standing controversy between an S = 1 model 
without an orbital degeneracy and the S =1/2 orbitally 
degenerate model of Castellani et al. [4]. In this Letter 
we propose a new model for the AF-ordered insulating 
(AFI) phase based on the molecular orbitals of the o 
axis V-V pairs, which combines features of both existing 
models and which reconciles the apparently conflicting 
experiments supporting each. 




FIG. 1. Corundum structure of V2O3. The V-ions 
(solid circles) are arranged in V-V pairs along the c-axis 
(face-sharing octahedra) and a honeycomb lattice in the other 
directions (edge-sharing octahedra). 

The V-ions in the corundum structure of V2O3 sit in a 
O-octahedron with a small trigonal distortion causing a 
small splitting in the non-bonding i2g-shell between the 
ai g -orbital oriented along the c-axis and doublet planar 
e ff -orbitals (see Fig. 1). In their early work, Castellani et 
al. proposed that one electron of the 3d 2 V 3+ -ion entered 
a spin singlet covalent ai 9 -bond in the V-V pair while 
the remaining electron was in the e g -doublet. Orbital 
ordering of these e s -doublets allowed them to explain 



the unusual magnetic structure of the AFI-phase with 
inequivalent n.n. exchange constants in the a — b plane 
(2 antiferromagnetic, 1 ferromagnetic) [5-7]. Paolasini et 
al. [8] interpreted their recent resonant x-ray experiments 
as a confirmation of this orbital ordering. On the other 
hand the polarized soft x-ray experiments by Park et al. 
[9] showed a coexistence of both (e g e g ) and (e g a\ g ) con- 
figurations in roughly equal amounts and these led Ezhov 
et al. [10] to argue for a S = 1 model with a (e g e g ) con- 
figuration and no orbital degeneracy. This is favored by 
the atomic Hund's Rule whose strength, as they point 
out, is not screened in the crystal. The differing planar 
exchange constants they attribute to the monoclinic dis- 
tortion in the AFI-phase. Yet general considerations of 
the phase diagram [11] and NMR investigations [12] all 
point towards to the presence of an orbital degeneracy. 

Here we take a different approach to the AFI-phase 
and start from an atomic limit but consider first the V-V 
pairs, since the intersite ai g -hopping matrix elements arc 
the largest [4]. Keeping a strong Hund's Rule coupling, 
as proposed by Ezhov et al., leads us to molecular orbitals 
for a V-V pair consisting of a superposition of [e g e g ) on 
one V-site and of (e g ai g ) on the second site with a total 
spin S = 2. This delocalized molecular orbital has also a 
two- fold degeneracy due to a choice in [e g a\ g ) among the 
e g -doublet. Next we consider planar hopping processes 
and show that in a reasonable parameter range the real 
spin (RS) structure is the most stable. This state has a 
ferro-arrangement of the molecular orbitals which agrees 
with the monoclinic structure and, as we shall sec, also 
with the x-ray experiments of Paolasini et al. 

Let us start with a description of a vertical pair. Fol- 
lowing Ref. [4], the two e g orbitals [13] are specified 
by a further index as |e ffl ) = l/y2(\d yz ) — \d zx )) and 
|e s a) = l/\fQ{2\d xy ) - \d yz ) - \d zx )), while the ai g orbital 

'xy) + \d yz ) + \d zx )). ForeachV 



is given by \ai g ) = 1/ V3(\d, 



ion, the d-orbitals are defined in a local coordinate sys- 
tem whose axis point towards the surrounding O ions, 
and thus refer to different, symmetry related orbitals for 
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the different V ions in the unit cell. Consequently, the e g \ 
and e g 2 orbitals on the two V ions of a vertical pair are 
not identical. This will be important when we compare 
our results to that of resonant scattering experiments. 
The intra-atomic interaction is described by three param- 
eters: U, the Coulomb interaction in the same orbital, U' , 
the Coulomb interaction in different orbitals, and J, the 
Hund's Rule coupling, which we assume satisfy the usual 
relation for t 2g orbitals: U = U' + 2 J. The trigonal crys- 
tal field induces an energy splitting A between the low- 
lying e g orbitals and the excited a\ g . Finally, the hopping 
integrals are denoted by i*- where 5 = a,b,c,d stands for 
the direction of the bond (a, b, c: bonds inside the hexag- 
onal planes, d: vertical bond) while — 1, 2, 3 denote 
the orbitals [e g \,e g2 and a\ g respectively). 

The main difference with Ref. [4] comes from the val- 
ues of the interaction parameters. The values used in 
Ref. [4] (U ~ 2eV, J = 0.2eV) are now believed to be 
much too small: Recent estimates based on LDA+U cal- 
culations [10] are in the range U ~ 5 cV and J ~ 1 eV. 
It turns out that this makes a dramatic difference for the 
ground state of a V-V pair. To be specific, if we con- 
sider the same hopping and crystal field parameters as in 
Ref. [4], and if we fix the ratio J/U = 0.1 to the value 
they used, there is a level crossing as a function of U 
between two very different situations. At small U, the 
ground state is 3-fold degenerate, with 3 levels nearby. 
This corresponds to the limit of Ref. [4] where two elec- 
trons go into the bonding molecular orbital built out of 
a\ g orbitals, the other two electrons being described by 
a spin 1/2 - pseudo spin 1/2 Kugel-Khomskii model [14]. 
At large U however, the ground state is 10-fold degen- 
erate. It corresponds to a total spin 2 with a two-fold 
degenerate orbital state. Since by symmetry tfj = if 
i 7^ j, this orbital wavefunction can actually be written 
down explicitly: 
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where e g stands for e g \ (e g2 ) in |— ) (|+)). This situation 
is generic for a large range of J/U including J/U = 0.2, 
and with the parameters proposed in Ref. [10], we found 
that the ground state is clearly of this second type. 

It is interesting to compare this state with the spin 1 
picture of Ezhov et al.. When the Hund's Rule coupling 
is large, all low-lying states can indeed be described by 
considering only the states with total spin 1 at each site. 
However, the resulting effective Hamiltonian for a pair of 
sites is not simply a Heisenberg Hamiltonian, since this 
would correspond to only 9 low-lying states. In fact there 
are 81 low-lying states, suggesting that it one wants to 
describe this system with a spin 1 operator, S, at each 
site, one should also include a pseudo-spin 1 operator, T, 
to describe the quasi-degeneracy of the t 2g orbitals. This 
orbital degree of freedom is crucial since it is responsible 
for a factor 2 in the 10-fold groundstate degeneracy. 

These results suggest that, instead of starting from a 
spin-orbital model with a spin 1/2 and a pseudo-spin 1/2 
at each V site, one should start from a spin-orbital model 



in which each vertical V-V pair is decribed by a spin 2 
for the total spin, say a, and a pseudo-spin 1/2 for the 
orbital degeneracy, say f, t z = 1/2 (r z = —1/2) corre- 
sponding to |+) (|— )) in Eq.l. The low energy properties 
are then determined by the way the degeneracy is lifted 
when these pairs are coupled by the in-plane hopping in- 
tegrals. Since these hopping parameters are small, we 
can treat them within second-order perturbation theory. 
For simplicity, we include only the largest hopping inte- 
gral t\\ = t, and the corresponding hopping integrals for 
directions b and c, in the present discussion. We have 
checked that the conclusions are unaffected by this sim- 
plification. The second order effective spin-orbital Hamil- 
tonian for n.n. along the a-axis then reads [15]: 



H cS (a) = Gai ■ a m + -G 3 (rf + t^)<t ; 



with 
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G = ~3 Gl + 3 G2 + 4 G3 ' Gl = A(U'-J) 



G 2 = 



4(U' + 2J) 



G, = 



t 2 



4(C/ + J) 



(3) 



The effective Hamiltonians for n.n. along the b- and 
c-axes are easily obtained by the trigonal rotation of 
H e ff (a) equivalent to the following replacement of the or- 
bital pseudo-spin t z -> -1/2t z ± V3/2t x . While there 
is a strong anisotropy in orbital space, the interaction 
preserves SU(2) symmetry for the spin operator a. 

Remarkably enough, the symmetry properties of this 
model are quite similar to the Kugel-Khomskii model for 
the cubic perovskite [14]. In fact one can regard the 
corundum lattice as a distorted simple cubic (sc) lattice 
of the V-V pairs. This analogy is useful to give a sys- 
tematic analysis for such a complicated system. Namely, 
it is promising that the stable magnetic phases of Kugcl- 
Khomskii-type models are collinear two-sublattice order- 
ings with associated orbital orderings. Within this cri- 
terion the possibility is naturally restricted into G, F, 
C, and A-type magnetic patterns [16]. In this language, 
the realistic magnetic structure of V2O3 corresponds to 
the G-type arrangement in the pair sc lattice: One of 3 
in-plane bonds is ferromagnetic and other two bonds are 
antiferromagnetic. 

Keeping in mind these relations, we have examined the 
stable phase in the molecular model by comparing the 
energies of all magnetic phases. This has been done, as 
for the Kugel-Khomskii model, within a mean-field de- 
coupling based on the order parameters (r"), (cr a ), and 
(r a cr' 3 ). Details will be given in a forthcoming paper 
[17]. The results are plotted in Fig. 2 as a function of 
Hund's Rule coupling J. It turns out that the stable 
phase changes successively from G to F phase as J in- 
creases. In order to gain energy by the orbital-dependent 
G3 term, the symmetry-broken C and A phases are sta- 
bilized in the intermediate- J region. In particular the 
realistic G-type phase is found to be the lowest for J/U 
around 0.2, which agrees with the estimates of Ezhov et 
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al., and which is consistent also with the stability region 
for the S = 2 degenerate molecular orbitals of a V-pair. 
For this phase, the orbital order parameter is ferromag- 
netic with t z < 0, i.e. the e g \ orbital is favoured (state 
| — ) of Eq. 1). Such a ferro molecular orbital order will 
cause an effective uniaxial stress on the lattice degrees of 
freedom, leading to a uniform rotation of V-V pairs. This 
is consistent with the monoclinic distortion proposed by 
Dernier and Marezio [20]. 
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FIG. 2. Energy comparison of various magnetic patterns. 
The letters refer to the notations of Ref. [16] for the sc lattice. 
The pictures give the corresponding magnetic pattern for the 
corundum lattice. The symbols |±), |+), and |— ) indicate 
which of the molecular orbitals of Eq. 1 is consistent with the 
magnetic pattern. The orbital degeneracy is lifted for states 
C and A, but not for states G and F. 

The physical picture that emerges from this model is 
very encouraging. First of all, the observed magnetic 
arrangement [5—7] is consistent with this model for rea- 
sonable values of the parameters. Second, the atomic 
configuration is a mixture of (e g e g ) and (ai s e g ), in agree- 
ment with X-ray absorption [9]. Third, there is an orbital 
degree of freedom whose ordering is consistent with the 
monoclinic distortion of the low-temperature phase [20]. 
It corresponds to choosing between e g \ and e g 2 for the 
V-V pairs. 

The results of our model are also consistent with the 
resonant x-ray scattering experiment of Paolasini et al. 
[8]. In that experiment, resonant scattering was observed 
in the low temperature phase at wavevector q = (111) 
and at energies corresponding to the transition from Is 
to unoccupied 3d states on V-ions. Since this is forbidden 
by symmetry if all the orbital configurations on V-atoms 
are equal, this led Paolasini et al. [8] to conclude to the 
existence of orbital ordering. They further suggested that 
their experiment can be taken to confirm one of the or- 
bital ordered phases previously proposed by Castellani et 
al [4]. We now show that the orbital ordering proposed 
in this paper for our model, while different from that in- 
terpreted by Paolasini et al., is also consistent with the 
experimental results. 

The resonant scattering amplitude as a function of the 



energy lu and q of x-rays for a crystal V2O3 is given by 

F = e ^' Pi fi(u), where fi(w) is the amplitude con- 

tributed from the V atom at position pi in the monoclinic 
unit cell of V2O3. The low-temperature monoclinic lat- 
tice of V2O3 has eight atoms in a unit cell (Fig. 3). Atoms 
1-4 have spin-up magnetic moments and atoms 5-8 have 
spin-down magnetic moments. Fat<j=(lll)is given by 

F111 - (/1 - h + h - U)e ia + (/ 2 - h + h - h)e~ ia 

where a is a phase factor which depends on p\ — P2- 
,fi{uj) depends in general on the magnetic moment and 
orbital occupation [18]. The nonvanishing intensity of 
(111) reflection for this energy implies that the combi- 
nations (/1 - / 5 + / 8 - fi) and (/ 2 - f 6 + f 7 - / 3 ) are 
nonzero. The fcrro-orbital phase in our model exhibits 
this feature for the following reasons. As discussed be- 
fore, the e g orbitals are defined with respect to a local 
co-ordinate system on each V ion. In particular, for the 
two V ions on a vertical bond, they are related by a ro- 
tation around the y-axis: C2(x,y, z) = (—x,y,—z) (the 
trigonal coordinate system is used here with z-axis di- 
rected perpendicular to the hexagonal plane) while for 
the V ions in the same hexagonal plane local coordinate 
systems are identical. Thus, the ferro-orbital phase actu- 
ally corresponds to having different orbitals on alternate 
hexagonal planes. Denoting these as 1 and 2 and denot- 
ing u and d for the spin up and down magnetic state, we 
then have e.g. 

fi - h + fa - U = /(!,«) - /(I, d) + /(2, d) - f(2, u) 

Since f depends on both the orbital occupation and the 
magnetic moment of the V-atom, j(u) ^ f(d). Thus, our 
model gives nonzero F m , and is qualitatively consistent 
with the experimental observation of Paolasini et al. [8] . 
More work is needed to compare our theory with the 
observed polarization and the azimuthal dependences of 
the resonance intensity. 

According to this explanation, the intensity of the 
(111) reflection is not simply a direct consequence of the 
orbital order, but comes both from magnetic and orbital 
order. This should be contrasted to Paolasini et al.'s 
explanation based on Ref. [4], where the form factor of 
Eq. (5) is non zero because the orbitals occupied on 
the two V of a vertical pair are different linear combi- 
nations of egi and eg2- While Castellani et al.'s picture 
is very specific to their S=l/2 model, there is room a 
priori within the S=l model for a similar orbital order- 
ing. Energetic considerations then show that the only 
serious candidate is an orbital ordering in which one V 
of a vertical pair would be in the (eg a\ g ) configuration, 
while the other one would be in the (eg eg) configura- 
tion. To be more precise, the model considered thus far 
corresponds to having all intcrplanc hopping amplitudes 
much larger than intraplane ones. If we maintain the 
limit of ig 3 large, but allows t\ x and t^ 2 to be compa- 
rable to intraplane hoppings, it can be shown that the 
problem can then be mapped into a transverse field Ising 
model. Details will be given in a forthcoming publica- 
tion [19]. Here we briefly summarize the main results. 
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In this mapping, there is an Ising spin associated with 
each V-V pair which corresponds to the orbital occupa- 
tion of (e g e g ) and (e g ai g ) on the two site. The transverse 
field strength h is given by the energy difference of the 
bonding and antibonding states of a V-V pair, a measure 
of the orbital quantum fluctuation in a V-V pair. The 
molecular orbital model corresponds to the large h limit. 
In the opposite limit h — ► 0, the spin and orbital order- 
ing depends on the relative strengths of the intraplane 
hoppings. In a reasonable parameter range, we found a 
ground state with RS spin and an orbital ordering corre- 
sponding to the pattern reported by Paolasini et al [8]. 
However, since in this state the V ions of a vertical pair 
are in (e g e g ) and (e g ai g ) configurations, the electronic 
densities are very different. This should lead to different 
local distortions of the O octahedra, which is inconsis- 
tent with the monoclinic structure reported by Dernier 
and Marezio [20], where all V are equivalent [21]. So we 
do not think that this kind of orbital ordering is realized 
in V 2 3 . 
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FIG. 3. The structure of the low temperature monoclinic 
phase of V2O3. The gray and filled circles correspond to 
spin-up and spin-down orientations of the local magnetic mo- 
ments on V ions. 

To summarize, we have proposed a model for the AF 
insulating phase of V2O3. This model seems to be the 
only way to combine basic facts about the electronic 
structure (S=l, orbital degeneracy, strong coupling along 
vertical pairs) into a coherent picture that agrees with all 
experiments. Further investigation of the resulting two- 
fold degenerate, S=2 model for the vertical pair is in 
progress. 
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